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Abstract 

We discuss a connection (and a proper place in this framework) of the un- 
forced and deterministically forced Burgers equation for local velocity fields 
of certain flows, with probabilistic solutions of the so-called Schrodinger inter- 
polation problem. The latter allows to reconstruct the microscopic dynamics 
of the system from the available probability density data, or the input-output 
statistics in the phenomenological situations. An issue of deducing the most 
likely dynamics (and matter transport) scenario from the given initial and ter- 
minal probability density data, appropriate e.g. for studying chaos in terms 
of densities, is here exemplified in conjunction with Born's statistical inter- 
pretation postulate in quantum theory, that yields stochastic processes which 
are compatible with the Schrodinger picture free quantum evolution. 



1 The Schrodinger reconstruction problem: most 
hkely microscopic dynamics from the input-output 
statistics data 

Probability measures, both invariant and nontrivially time-dependent, often on dif- 
ferent levels of abstraction, are ubiquituous in diverse areas of physics. According 
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to pedestrian intuitions, one normally expects that any kind of time develope- 
ment (dynamics, be it deterministic or random), which is analyzable in terms of 
probability, under suitable mathematical restrictions may give rise to a well defined 
stochastic process. Non-Markovian implementations are regarded to be close to real- 
ity, but the corresponding Markovian approximations (when appropriate) are easier 
to handle analytically. 

Given a dynamical law of motion (for a particle as example), in many cases 
one can associate with it (compute or approximate the observed frequency data) 
a probability distribution and various mean values. In fact, it is well known that 
inequivalent finite difference random motion problems may give rise to the same 
continuous approximant (like e.g. in case of the diffusion equation representation of 
discrete processes). As well, in the study of nonlinear dynamical systems, specifi- 
cally those exhibiting the so-called deterministic chaos, 0, H, |], given almost any 
(basically one- dimensional in the cited references) probability density, it is possible 
to construct an infinite number of deterministic finite difference equations, whose 
iterates are chaotic and which give rise to this a priori prescribed density. 
Studying dynamics in terms of densities of probability measures instead of individ- 
ual paths (trajectories) of a physical system is a respectable tool, |^, even if we 
know exactly the pertinent microscopic dynamics. 

Under general circumstances, the main task of a physicist is to fit a concrete dynami- 
cal model (through a clever guess or else) to available phenomenological data. Then, 
the distinction between the chaotic (nonlinear, deterministic) and purely stochas- 
tic implementations may not be sharp enough to allow for a clean discrimination 
between those options: the intrinsic interplay between the stochastic and determin- 
istic modelling of physical phenomena, 0, blurs the access to reality and certainly 
precludes a definitive choice of one type of modelling against the other. 

An inverse operation of deducing the detailed (possibly individual, microscopic) 
dynamics, which is either compatible with a given probability measure (we shall be 
mostly interested in those admitting densities) or induces its own time evolution, 
cannot have a unique solution. However, the level of ambiguities can be substan- 
tially reduced, if we invoke the so-called Schrodinger problem of reconstructing the 
microscopic dynamics from the given input-output statistics data and/or from the 
a priori known time developement of a given probability density. The problem is 
known to give rise to a particular class of solutions {most likely interpolations), in 
terms of Markov diffusion processes, 0-0. 

In its original formulation, due to Schrodinger, 1^-0, one seeks the answer to the 
following question: 

given two strictly positive (usually on an open space-interval) boundary probability 
densities Pq{x), Pt{x) for a process with the time of duration T > 0. Can we uniquely 
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identify the stochastic process interpolating between them ? 

Another version of the same problem, @ , departs from a given (Fokker- Planck- 
type) probability density evolution and investigates the circumstances allowing to 
deduce a unique random process from this dynamics. We shall pay some attention 
to this issue in Section 3. 

The answer to the above Schrodinger's quaestion is known to be affirmative, if 
we assume the interpolating process to be Markovian. In particular, we can get 
here a unique Markovian diffusion process which is specified by the joint probability 
distribution 

mT{A,B)= / d^x I d^ymT{x,y} (1) 

J A J B 

(fymT{x,y) = po{x) 

(fxmT{x,y) = priy) 
where 

rnrix, y) = Uo{x) k{x, 0, y, T) vriy) (2) 

and the two unknown functions Uq{x)^ VT{y) come out as solutions of the same sign 
of the integral identities (1). Provided, we have at our disposal a continuous bounded 
strictly positive (ways to relax this assumption were discusssed in Ref. [Q) integral 
kernel s,y,t),0 < s < t < T. 

We shall confine further attention to cases governed by the familiar Feynman- 
Kac kernels. Then, the solution of the Schrodinger boundary-data problem in terms 
of the interpolating Markovian diffusion process is found to be completely specified 
by the adjoint pair of parabolic equations. In case of gradient forward drift fields, 
the pertinent process can be determined by checking (this imposes limitations on 
the admissible potential) whether the Feynman-Kac kernel 

k{y, s, x,t) = J exp[- c(c^(r), T)dT]dp^(^''J^{uj) (3) 

is positive and continuous in the open space-time area of interest (then, additional 
limitations on the path measure need to be introduced, |11[), and whether it gives 
rise to positive solutions of the adjoint pair of generalised heat equations: 

dtu{x, t) = h'Au{x, t) — c(x, t)u{x, t) (4) 

dtv{x, t) = —uAv{x, t) + c(x, t)v{x, t) . 

Here, a function c{x, t) is restricted only by the positivity and continuity demand 
for the kernel (3), see e.g. |jlO|- In the above, dp%^}{uii) is the conditional Wiener 



measure over sample paths of the standard Brownian motion. 
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Solutions of (4), upon suitable normalisation give rise to the Markovian diffusion 
process with the factorised probability density p{x,t) = u{x,t)v{x,t) which, while 
evolving in time, interpolates between the boundary density data p{x, 0) and p{x, T). 
The interpolation admits a realisation in terms of Markovian diffusion processes with 
the respective forward and backward drifts defined as follows: 

b{x,t)=2u (5) 

V{X, t) 

K{x,t) = -2u 

U[X, t) 

in the prescribed time interval [0,T]. 

The related transport equations for the densities easily follow. For the forward 
interpolation, the familiar Fokker-Planck equation holds true: 

dtp{x, t) = z/Ap(x, t) — V[6(x, t)p(x, t)] (6) 

while for the backward interpolation we have: 

dtp{x,t) = -uAp{x,t) - V[K{x,t)p{x,t)] . (7) 

We have assumed that drifts are gradient fields, curlh = 0. As a consequence, 
those that are allowed by the prescribed choice of c(x, t) must fulfill the compatibility 
condition 

c(f,t) = 9,$ + ^(|; + V6) (8) 

which establishes the Girsanov-type connection of the forward drift t) = 2i'V^{x, t) 
with the Feynman-Kac, cf. |]1T], |T0|, potential c{x, t). In the considered Schrodinger's 
interpolation framework, the forward and backward drift fields are connected by the 
identity 6* = 6 — 2uVlnp. 

One of the distinctive features of Markovian diffusion processes with the posi- 
tive density p{x,t) is that, given the transition probability density of the (forward) 
process, the notion of the backward transition probability density p^{y, s, x,t) can 
be consistently introduced on each finite time interval, say < s < t < T: 

p{x, t)p^{y, s, X, t) = p{y, s, f, t)p{y, s) (9) 

so that / p{y, s)p{y, s, x, t)d^y = p{x, t) and p{y, s) = J p^{y, s, x, t)p{x, t)d^x. 

The transport (density evolution) equations (6) and (7) refer to processes running 
in opposite directions in a fixed, common for both, time-duration period. The 
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forward one, (6), executes an interpolation from the Borel set A to B, while the 
backward one, (7), executes an interpolation from B to A, compare e.g. the defining 
identities (1). 

The knowledge of the Feynman-Kac kernel (3) implies that the transition prob- 
ability density of the forward process reads: 

p{y, g, f , t) = k{y, s, X, t) ^^Z' ^\ ■ (10) 

v{y, s) 

while the corresponding (derivable from (10), since p{x,t) is given) transition prob- 
ability density of the backward process has the form: 

u(y^ ts) 

P*{y, s, X, t) = k{y, s, X, t) ' . (11) 

u{x^ t) 



Obviously, [T^, |^, in the time interval < s < t < T there holds: 

u{x,t) = / uo{y)k{y,s,x,t)d^y (12) 



v{y, s) = J k{y, s, x, T)vT{x)d^x . 

Consequently, the system (4) fully determines the underlying random motions, for- 
ward and backward, respectively. 



2 The Burgers equation in Schrodinger's inter- 
polation 

The prototype nonlinear field equation named the Burgers or "nonlinear diffusion" 



equation (typically without, fT^, |T5[, the forcing term F{x,t)): 



dtVB + ivBV)vB = lyAvB + F{x, t) (13) 
recently has acquired a considerable popularity in the variety of physical contexts. 

By dropping the force term in (13), we are left with a commonly used form of the 
"nonlinear diffusion equation" whose solutions are known exactly, in view of the 
Hopf-Cole linearising transformation mapping (13) into the heat equation. Here, 
dtVB + {vbV)vb = uAvb is mapped into dt9 = uAO, by means of the substitution 
vb = —2i>Vln6. This linearisation of the Burgers equation is normally regarded to 
be devoid of any deeper physical meaning, and specificaly the link with stochastic 



5 



processes determined by the heat equation has not received a proper attention. Our 
previous analysis shows that the intrinsic interplay between the deterministic and 
random evolution, appropriate for a large class of classically chaotic systems, extends 
to much wider framework. 

Burgers velocity fields can be analysed on their own with different (including 
random) choices of the initial data and/or force fields. However, we are interested 
in the possible diffusive matter transport that is locally governed by Burgers flows, 
cf. [|13[. In this particular connection, let us point out a conspicouous hesitation 



that could have been observed in attempts to establish the most appropriate matter 
transport rule, if any diffusion-type microscopic dynamics assumption is adopted to 
underly the "nonlinear diffusion" (13). 

Depending on the particular phenomeno logical departure point, one either adopts 



the standard continuity equation, |Tq, 17|, that is certainly valid to a high degree of 



accuracy in the so-called low viscosity limit z/ | , but incorrect on mathematical 
grounds if there is a genuine Markovian diffusion process involved and simulta- 
neously a solution of (13) stands for the respective current velocity of the flow: 
dtp{x,t) = -V[v{x,t)p{x,t)] . 

Alternatively, following the white noise calculus tradition telling that the stochas- 
tic integral X(t) = Jq vb{X{s), s)ds + /q f]{s)ds necessarily implies the Fokker-Planck 
equation, one is tempted to adopt: dtp{x,t) = z/Ap(x, t) — V[vB{x,t)p{x,t)] which 
is clearly problematic in view of the classic Mc Kean's discussion of the propagation 



of chaos for the Burgers equation, ITq, 19, 20| and the derivation of the stochastic 



"Burgers process" in this context: "the fun begins in trying to describe this Burgers 
motion as the path of a tagged molecule in an infinite bath of like molecules" 



To put things on the solid ground, let us consider a Markovian diffusion process, 
which is characterised by the transition probability density (generally inhomoge- 
neous in space and time law of random displacements) p{y, s,x,t) , < s < t < T, 
and the probability density p{x,t) of its random variable X{t) , < t < T. The 
process is completely determined by these data. For clarity of discussion, we do not 
impose any spatial boundary restrictions, nor fix any concrete limiting value of T 
which, in principle, can be moved to infinity. 

Let us confine attention to processes defined by the standard backward diffu- 
sion equation. Under suitable restrictions (boundedness of involved functions, their 
continuous differentiability) the function: 



g{x, s) = J p{x, s, y, T)g{y, T)d y (14) 

satisfies the equation 

-dsg{x, s) = i^Ag^x, s) + [b{x, s)V]g{x, s) . (15) 
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Let us point out that the vahdity of (14) is known to be a necessary condition for 
the existence of a Markov diffusion process, whose probabihty density t) is to 
obey the Fokker-Planck equation (the forward drift 6(x, t) replaces the previously 
utilized Burgers velocity VB{x,t)). 

The case of particular interest, in the traditional nonequilibrium statistical physics 
literature, appears when p{y,s,x,t) is a fundamental solution of (15) with respect 
to variables y,s, [^, see however |TD| for an analysis of alternative situa- 

tions. Then, the transition probability density satisfies also the second Kolmogorov 
(e.g. the Fokker-Planck) equation in the remaining x, t pair of variables. Let us 
emphasize that these two equations form an adjoint pair of partial differential equa- 
tions, referring to the slightly counterintuitive for physicists, though transparent for 
mathematicians, ^ ^ , issue of time reversal of diffusions. 

We can consistently introduce the random variable of the process in the form 
X{t) = /q b{X{s), s) ds + \/2h'W{t). Then, in view of the standard rules of the Ito 
stochastic calculus, ||28|, ^ P|, we realise that for any smooth function /(x, t) of the 
random variable X{t) the conditional expectation value: 

lim^tio^^ij p{x,t,y,t + /\t)f{y,t + /\t)d''y- f{x,t)] = {D+f){X{t),t) = (16) 

= {dt + hV + v/\)f{x,t) , 

where X{t) = x, determines the forward drift b{x,t) of the process (if we set com- 
ponents of X instead of /) and, moreover, allows to introduce the local field of 
(forward) accelerations associated with the diffusion process, which we constrain by 
demanding (see e.g. Refs. [^, |2^, ^ for prototypes of such dynamical constraints): 

{DlX){t) = {Dj){X{t),t) = {dtb+ (bV)b + uAb){x,t) = F{x,t) (17) 

where X{t) = x and, at the moment arbitrary, function F{x,t) may be interpreted 
as an external forcing applied to the diffusing system, |[TT[] . 

By invoking (9), we can also define the backward derivative of the process in the 



conditional mean (cf. |]TT], |2^, ^ for a discussion of these concepts in case of the 



most traditional Brownian motion and Smoluchowski-type diffusion processes) 

l^mAtlo ^[^- / PMt- At,x,t)0y] = {D^X){t) = K{X{t),t) (18) 

(/}_/) (X(t),t) = {dt + IV - uA)f{X(t),t) 
Accordingly, the backward version of the acceleration field reads 

{DlX){t) = {Dlxm = F{X{t\t) (19) 
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where in view of 6=,, = b — 2v'S/lnp we have exphcitly fulfilled the forced Burgers 
equation : 

dtK + - v^K = F (20) 

and, |]2^, H, |TT|, under the gradient-drift field assumption, curlh^ = 0, we deal with 
F{x, t) = 2z/Vc(x, t) where the Feynman-Kac potential (3) is explicitly involved. 

Let us notice that the familiar (linearization of the nonlinear problem) Hopf- 
Cole transformation, [|15], 0, of the Burgers equation into the generalised diffusion 
equation (yielding explicit solutions in the unforced case) has been explicitly used 
before (the second formula (4)) in the framework of the Schrodinger interpolation 
problem. In fact, by defining = logu, we immediately recover the traditional 
form of the Hopf-Cole transformation for Burgers velocity fields: 6* = — 2z/V$*. 
In the standard considerations that allows to map a nonlinear (unforced Burgers) 
equation into a linear, heat, equation. In the special case of the standard free 
Brownian motion, there holds b{x, t) = while 6=k(x, t) = —2v'Wlog p{x, t). 

Let us point out that the equation (7) is in fact the only transport equation 
where the Burgers velocity field is allowed to be undisputably present, under the 



diffusive scenario assumption, |T3[. The standard continuity equation is certainly 



inappropriate for nonzero values od the diffusion constant u. 



3 Reconstruction of the microscopic dynamics from 
the probabihty density data: obstacles exempli- 
fied 

We have mentioned before, that another version of the Schrodinger boundary data 
problem, , departs directly from a given (Fokker-Planck-type) probability density 
evolution and investigates the circumstances allowing to deduce a unique random 
process from this dynamics. Surely, solutions of the Fokker-Plack equation itself do 
not yet determine the underlying stochastic process. Additional assumptions are 
always necessary and a number of traps to be avoided. 

As a particular guide to those obstacles, we shall refer to the familiar free quan- 
tum evolution that is regarded as the time adjoint parabolic problem, exactly in the 
spirit of our previous discussion. 

In our previous paper, the major conclusion was that in order to give a 



definitive probabilistic description of the quantum dynamics as a unique diffusion 
process solving Schrodinger's interpolation problem, a suitable Feynman-Kac semi- 
group must be singled out. Let us point out that the measure preserving dynamics. 
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permitted in the presence of conservative force fields, was investigated in 111 



see 



also f?, 32 



The present analysis was performed quite generally and extends to the dynam- 
ics affected by time dependent external potentials, with no clear-cut discrimination 
between the nonequilibrium statistical physics and essentially quantum evolutions. 
The formalism of Section 1 encompasses both groups of problems. Nevertheless, it 
is quite illuminating to see directly how sensitive, even in simplest cases, the for- 
malism is with respect to any attempt of relaxing our previous assumptions and the 
Schrodinger interpolation problem rules-of-the-game. Specifically in the quantum 
domain, where the seemingly trivial case of the free evolution, which is nonstation- 
ary, needs the general parabolic system (4) to be considered. Even worse, then the 
system (4) displays a nontrivial nonlinearity: the parabolic equations are coupled 
by the effective, solution dependent potential. At the first glance, this feature might 
seem to exclude the existence of any conceivable Feynman-Kac (dynamical semi- 
group) kernel, and in consequence any common-sense law of random displacements 
(i.e. the transition probability density) governing the pertinent stochastic evolu- 
tion. Certainly, the existence of fundamental solutions in this case is far from being 
obvious. 

At this point, let us emphasize that our principal goal is to take seriously the 
Schrodinger picture quantum dynamics under the premises of the Born statistical 
postulate. Hence, once we select as appropriate a concrete quantal interpolation 
between the prescribed (phenomenologically supported in particular) input-output 
statistics data po(^) and Pt{x) in terms of p{x,t) = 4'{x,t)ip{x,t), t G [0,T], where 
ip{x,t) solves the Schrodinger equation then, on exactly the same footing, we are 
entitled to look for an alternative probabilistic explanation (or appropriate descrip- 
tion) of the very same interpolation, in terms of a well defined Markov stochastic 
(eventually diffusion) process. 

We shall proceed in the spirit of Section 1, while restricting our discussion to 



the free Schrodinger dynamics. Following Ref. |^ we shall discuss the rescaled 
problem so as to eliminate all dimensional constants. 

The free Schrodinger evolution idfip = —Aip implies the following propagation 
of a specific Gaussian wave packet: 

^(x,0) = (2vr)-V4exj9(-^) ^ (21) 

^{x,t) = {-fl' {2 + 2^tr^''exp[--^^] 
vr 4(^1 + it) 

So that ^ 

Po{x) = \^{xM^ = (27r)"i/2gxp[-y] (22) 



9 



p{x,t)^\4^{x,t)f^[27ril + f exp[- 



2(l+t2)J 

and the Fokker-Planck equation (easily derivable from the standard continuity equa- 
tion dtp — — V(vp), v{x,t) = xt/{l + t^)) holds true: 

dtp^Ap-Vibp) , b{x,t)^-^^x (23) 

The Madelung factorization ip = exp{R + iS) implies (notice that v = 2VS 
and b = 2V(-R + S)) that the related real functions 9{x,t) = exp{R + S) and 
9^{x, t) = exp{R — S) read: 

9{x,t) = [27r(l + V4ea;p(-^ - -arctant) 

9,ix,t) = [27vil + t^)]-'/'exp{-^^^ + larctant) (24) 
They solve a suitable version of the general parabolic equations (4), namely : 

dt9 = -A9 + (25) 
dt9^ = A9i, — —fl9^ 

with 

1 x^ \ A 0^/2 

-Q(x, t) = , — = 2^=f— = Q(x, t) (26) 

2 ^ ' ^ 2(l + t2)2 i+t2 pi/2 > y V y 

By setting t = T we associate with the above dynamics the terminal density 
Pt{x), and then the concrete Schrodinger boundary data problem for the stochastic 
interpolation po{x) — > Pt{x), (1). 

To capture the spirit of our previous discussion, wc shall replace equations (25) 
by the more general equations (4), where only the potential c{x,t) will be identified 
with the above jfl{x,t). Then, we shall look for solutions u{x,t), v{x,t) of these 
parabolic equations, and in particular we shall identify the quantally implemented 
functions ^*(x, t), 9{x, t), (24), amon^' them. Effectively, it amounts to the previously 
mentioned linearisation of the nonlinear parabolic system. 

In view of the relatively simple form of the probability density p{x,t), (22) one 
might be tempted to guess (more or less fortunately) the transition probability 
density, consistent with the propagation (22). However, it is well known that there 
are many stochastic processes implying (22) for all t e [0, T], which not necessarily 
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have much in common with the original wave function dynamics (21). In general they 
are incompatible with the corresponding parabolic system (cf. (4) and (25)). If it 
happens otherwise, the reason for this proliferation of would-be consistent stochastic 
processes is rooted in exploting the particular functional form of solutions, instead 
of relying on the form-independent arguments, e.g. (4). 

Let us consider simple examples which, albeit coming under very special cir- 
cumstances (free dynamics with a specific initial wave packet choice, and no zeros 
admitted in the course of the propagation), clearly idicate how important is the 
proper choice of the Feynman-Kac kernel. The virtue of a parabolic system (4) is 
that its form is universal for the Schrodinger dynamics, and thus does not depend on 
a particular functional form of solutions nor this of external potentials. It appears 
that the system (4) sets a very rigid framework for the probabilistic manifestations 
(e.g. stochastic processes) of the quantum Schrodinger dynamics. 



Example 1: We shall demonstrate that an improper {not through (4) or (25)), 
but fortunate, choice of the kernel might lead to an alternative stochastic represen- 
tation of the quantum dynamics (22). 

Let us begin from directly introducing the transition probability density 

p{y, s, X, t) = [2n{t' - s')]-'/' exp[-^f^] (27) 

which for all intermediate times < s < t < T executes a desired propagation 
p{x,t) — J p{y,s,x,t)p{y,s)dy, (22). Clearly, the Chapman-Kolmogorov identity 
/ p{y, s, z, t)p{z, t, X, t)dT — p{y, s, x, t) holds true, and the properties (the first one, 
for all e > 0): 

l^IPrv TTI / Piv^ t,x,t + At)dx = 

AUO At J\x-y\>e 
1 r+oo 

limAtio-jTL {x-y)p{y,t,x,t + At)dx = (28) 

ZAt J —oo 

lirriAm-r- I {x - y) p{y,t,x,t + At)dx = 2t 

ZAt J — oo 

tell us that the law of random displacements p{y, s, x, t), (27), can be attributed to 
a Markov diffusion process associated with the parabolic (Fokker- Planck) equation 

dtp = tA^p (29) 

In fact, our p{y, s, x, t) is a fundamental solution of this equation with respect to x, t 
variables, while obeying the time adjoint parabolic equation in the remaining (e.g. 
y, s) pair of variables 

dsp{y, s, X, t) = -sAyp{y, s, x, t) (30) 
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This diffusion has a vanishing forward drift and the quadratic in time variance 
(the diffusion coefficient equals t), hence its local characteristics are completely 
divorced from those of the Nelson process [PD| derivable from the solution (21) of 
the Schrodinger equation. 

Interestingly, since p{y,s,x,t) itself is a perfect, strictly positive and continuous 
in all variables (Markov) semigroup kernel, nothing prevents us from performing 
the Schrodinger problem analysis (1) with the boundary densities po{x) and Pt{x) 
defined by the above free evolution problem. However, we shall proceed otherwise 
and having given explicit solutions of the parabolic system (25) we introduce another 
strictly positive and continuous in all variables function: 



h{y,s,x,t)=p{y,s,x,tf-^ = [27r(t2 - .^)]-V2(l±l.)i/4 (31) 



. (x - y)^ . . y"^ 1 - s x'^l-t. A. 

Gxp\ — —\ expl H expl-iarctan t — arctan s)\ 

2(t2-s2)J 4 1 + s2 4 1 + t2J ^^2^ 

and observe that the Schrodinger system (1) in the present situation is involved as 
well, since trivially there holds: 

Po{x) = e^x, 0) J h{x, 0, y, T)e{y, T)dy (32) 

PT{x) = e{x,T) J hiy,0,x,T)9.{y,0)dy 

Disregarding the derivation which has led us to (22), we can simply consider (22) 
as the Schrodinger system of equations with a fixed kernel and boundary density 
data. Then, we immediately infer that by Jamison's theorem, 0, its unique (up 
to a coefficient) solution is constituted by the pair 9^{x,0), 9{x,T) of functions, 
already determined by (24). Moreover, ki{y, s, x, t) obeys the Chapman-Kolmogorov 
composition rule: 

ki{y, z,T)ki{z,T^x, .t)dT = (33) 



In view of / p{y, s, x, t)dx = 1 for all s < t, we have 

J k,{x,s,y,t)e{y,t)dy = e{x,s) (34) 

and, since ^6** = p, we get 

/0(,y, s) 
0*{y,s)p{y,s,x,t)jj-^dy = (35) 



12 



Thus, undoubtedly we have in hands a complete solution of the Schrodinger bound- 
ary data problem (1): for the once chosen kernel fci, this solution is unique, and com- 
patible with the dynamics of the corresponding Schrodinger wave function. But, the 
constructed stochastic process is completely incongreuent with the standard wisdom 
about Nelson's diffusion processes P, ^ |30|. The reason is clear: our analysis 
was performed for a particular solution, whose functional form allows for an alter- 
native stochastic representation. But, let us stress the point, if we look for the 
functional-form-independent construction, it is the parabolic system (4) from which 
one should depart. 

Anyway, even the inappropriate choice of the integral kernel ki , does allow to de- 
rive the quantum mechanically implemented dynamics (22) from respectively 9{x, T) 
and 6'*(x, 0), by means of the propagation formulas (4). The probability density 
evolves in time correctly, but the vanishing drift and the linear in time diffusion 
coefficient situate this stochastic process outside the scope set by (25) and (4). 



Example 2: We shall demonstrate, that another choice of the kernel, still with 
no reference to the system (4), will allow to reproduce the stochastic propagation 
with the probability density, drifts and diffusion coefficient of Nelson's stochastic 
mechanics, which however is not Nelson's process for the quantum evolution (22). 



We are inspired by our previous paper |30], where an interesting stochastic propa- 



gation, compatible with (22), was introduced by means of the transition probabihty 
density: 

p,,.(x,t) = [47r(t - s)]-V2 exM-^^^^^] 

Py,s{x, s) = 6{x — y) , < s < t < T (36) 



^ 1 + S2 ' 

Here, the density p{y, s), (22), is propagated into the corresponding p(x, t) according 
to the rule p{x,t) = J Py^s{x,t)p{y, s)dy, for all intermediate times < s < t < T. 
As noticed in [0, this propagation is somewhat pathological since it does not obey 
the Chapman-Kolmogorov composition rule: / Py^s{,z,T)pz,T{,x,t)dT ^ py^s{x,t) and 
thus p cannot be interpreted as a transition density of the Markov process. 

However, if we would naively proceed like in the Example 1 and define the strictly 
positive continuous function 

hiy, s, X, t) = py^s{x, t)^f^ (37) 



13 



where < s < t < T and 6{x,t) is given by (24), then the Schrodinger system 
(32), with ^2 replacing ki, trivially appears. Indeed, because of / Py^s{x,t)dx = 1 
for s < t, there holds: 



e,{x,0) I k2{x,0,y,T)e{y,T)dy = e,{x,0) J p,^o{y,T)e{x,0)dy = 

= e,{x,o)eix,o) = po{x) 

e{x, T) I k2{y, 0, X, T)9,{y, 0)dy = f Py^ix, T)e{y, 0)^,(t/, Q)dy = (3^ 



Pyfi{x, T)p{y, 0)dy = pt{x) 

As a consequence, if we analyze the above Schrodinger system with the boundary 
data Pq{x) and Pt{,x) fixed by (22) (as before), but with the new kernel /c2 then, 
somewhat unexpectedly, the same as before pair 9{x,Q),9^{x,T) necessarily comes 
out as a solution. Let us emphasize that the solution is unique for the chosen kernel 
/c2, albeit it coincides with the unique (as well) solution previously associated with 
the kernel ki (cf. Example 1). 

The meaning of the uniqueness of solution of the Schrodinger system ^ becomes 
clear: if we have prescribed the boundary density data the solution is unique for 
a chosen kernel, but there are many kernels which may give rise to the very same 
solution. 

The pathology (non-Markovian density) of Py^s{x,t) extends to k2{y, s,x,t) and 
the semigroup composition rule is invalid in this case. Nevertheless, we can blindly 
repeat the step (32), with k2 instead of ki, so reproducing the evolution (22). More- 
over, in the present case, [0, we can exploit the standard recipe to evaluate the 
forward drift of a conventional diffusion: 

If 1 — t 

limAtio^lJ yPxAv^t + ^t)dy - x] = b{x,t) = - ^^^^ x (39) 



Clearly, it is the forward drift of the Nelson diffusion p7| , |30| associated with (24) 



and it consistently appears in the corresponding Fokker-Planck equation (6). 

Let us observe that Py,s{x,t) solves the first Kolmogorov equation with respect 
to X, t: 

dtPy,s{x, t) = A^Py^s{x, t) - by^s(t)VxPy,s{^, t) (40) 



byA't) = y 



dct,s 
dt 
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As such, it can be exploited to construct a genuine Markov process, albeit discon- 
nected from the quantal dynamics (22). Namely, we can define another solution of 
the equation (40), in variables 

Py,sixi,ti,X2,t2) = [A7T{t2-h)]-'/^exp[-^^^—^^^-^] (41) 

C = Ct2,s - , 0<S<ti<t2<T 

with Ct^s given by (36). It is easy to verify that the transition density (41) actually 
is a fundamental solution, and as such satisfies the second Kolmogorov equation 
with respect to X2,t2, for each fixed y,s label, 0<s<ti<t2<T. Consequently, 
we have in hands the {y, s)-family of well defined Markovian transition probability 
densities Py^s for random propagation scenarios. Indeed, to this end one needs to 
check the (apparent) compatibility conditions: (a) Py,s{xi,t,X2,t) = 6{x2 — xi), 

(b) J Py,s{xi,ti,X2,t2)py,s{xi,ti)dxi =Py,s{x2,t2) and in addition 

(c) J Py,s{xi,ti,X2,t2)py,s{x2,t2,X3,t3)dx2 = Py^s{xi,ti, X2.,h) , wherc < s < ti < 
h < ts < T and Py^s{x,t) plays the role of the density of the Markov process. The 
identity (c) in the above is the Chapman- Kolmogorov formula. 

To avoid the above obstacles, the only how-to-proceed procedure is provided by 
the route outlined before, e.g. that leading from the Feynman-Kac kernel to the 
associated Markov diffusion process via Schrodinger's boundary-data problem. The 
complete solution to this particular issue, in the quantum dynamics context, has 
been given elswhere, . 
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